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ABSTRACT 

We present densely sampled BVRI light curves of the optical transient associated with the gamma-ray burst 
GRB 030329, the result of a coordinated observing campaign conducted at five observatories. Augmented with 
published observations of this GRB, the compiled optical dataset contains 2687 photometric measurements, 
obtained between 78 minutes and 79 days after the burst. This dataset allows us to follow the photometric 
evolution of the transient with unprecedented detail. We use the data to constrain the light curve of the under- 
lying supernova 2003dh, and show that it evolved faster than, and was probably somewhat fainter than the type 
Ic SN 1998bw, associated with GRB 980425. We find that our data can be described by a broken power-law 
decay perturbed by a complex variable component. The early- and late-time decay slopes are determined to 
be ai « LI and a2 ~ 2. Assuming this single-break power-law model, we constrain the break to lie between 
3 and ^ 8 days after the burst. This simple, singly-broken power-law model, derived only from the analysis 
of our optical observations, may also account for available multi-band data, provided that the break happened 
^ 8 days after the burst. The more complex double-jet model of Berger et al. provides a comparable fit to the 
optical. X-ray, mm and radio observations of this event. The unique early coverage available for this event al- 
lows us to trace the color evolution of the afterglow during the first hours after the burst. We detect a significant 
change in optical colors during the first day. Our color analysis is consistent with a cooling break frequency 
sweeping through the optical band during the first day. The light curves of GRB 030329 reveal a rich array 
of variations, superposed over the mean power-law decay. We find that the early variations (< 8 days after 
the burst) are asymmetric, with a steep rise followed by a relatively slower (by a factor of about two) decline. 
The variations maintain a similar time scale during the first four days, and then get significantly longer The 
structure of these variations is similar to those previously detected in the afterglows of several GRBs. 

Subject headings: Gamma Rays: Bursts — Stars: Supernovae: General — Stars: Supernovae: Individual: 
Alphanumeric: SN 2003dh 



I. introduction 

Following the discovery of low-energy transients associ- 
ated with long duration (> 2 s) gamma-ray bursts (GRBs; van 
Paradijs et al. 1997), a major effort was made to characterize 
the temporal evolution of these sources across the electromag- 
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netic spectrum. In the optical regime, the associated transient 
sources were found to decline rapidly with time. The emis- 
sion from the first optical transients (OTs) discovered was re- 
ported to decay as a power-law with time, extending from the 
epoch of discovery and continuing for tens of days (e.g., GRB 
970228 - Galama et al. 1997; GRB 970508 - Galama et al. 
1998a; Sokolov et al. 1998; Bloom et al. 1998). 

The decay slopes measured for OTs were used to constrain 
explosion models. In particular, with the increasing popular- 
ity of non-isotropic models (involving highly relativistic jets 
or cones), a temporal break in the decline slope was predicted 
(e.g., Rhoads 1997; Panaitescu, Meszaros, & Rees 1998; Sari, 
Piran & Halpern 1999) - and detected (e.g., Castro-Tirado et 
al. 1999; Kulkai-ni et al. 1999; Stanek et al. 1999; HaiTison et 
al. 1999; Price et al. 2001). 

Following the discovery of supernova (SN) 1998bw in the 
eiTor box of GRB 980425 (Galama et al. 1998b) the associa- 
tion of GRBs with SN explosions came into focus. Late-time 
"bumps" in OT light curves were interpreted as the signature 
of underlying SN explosions (e.g.. Bloom et al. 1999). While 
observational evidence supporting the SN hypothesis accumu- 
lated (e.g.. Bloom et al. 2002; Garnavich et al. 2003; Price 
et. al. 2003a), a direct observational proof for the existence of 
an underlying SN explosion has long remained missing, and 
alternative explanations for the origin of these bumps were 
suggested (e.g., Waxman & Draine 2000; Esin & Blandford 
2000; Reichart 2001). 

With growing interest and efforts by the astronomical com- 
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munity, the observed OT light curves became increasingly 
better-measured. In particular, for some OTs, a dense tem- 
poral sampling of the light curve, sometimes starting shortly 
(minutes to hours) after the GRB trigger, was carried out by 
world-wide observing networks. A case to point is the OT of 
GRB 021004, where multiple bumps and wiggles in the light 
curve, probably not associated with a SN, were first observed 
(e.g., Bersier et al. 2003; Mirabal et al. 2003; Fox et al. 2003). 

On 2003 March 29, a bright GRB was detected by the 
HETE-II spacecraft (Vanderspek et al. 2003). The early dis- 
covery of the associated OT (Peterson & Price 2003; Torii 
2003; Torii et al. 2003; Price et al. 2003b; Uemura et al. 
2003; Sato et al. 2003), and its brightness, triggered a world- 
wide observational effort, involving tens of ground and space- 
based facilities, observing in various wavelengths, from X-ray 
to radio. The brightness of the OT allowed the prompt deter- 
mination of the redshift of the source. At z=0.1685 (Greiner 
et al. 2003a) this event is the closest GRB to date for which 
a typical OT was discovered. The relatively low redshift 
of this burst presented a unique opportunity to search for a 
clear spectroscopic signature of an underlying SN. Indeed, 
intensive spectroscopic monitoring of the optical source re- 
vealed the emerging spectrum of a type Ic SN 1998bw-like 
event, designated SN 2003dh, conclusively proving that at 
least some of the long-duration GRBs are associated with SN 
explosions (Stanek et al. 2003; Hjorth et al. 2003; Matheson 
et al. 2003). 

Numerous works have akeady presented observations of 
this unique event. Early optical observations, obtained with 
small telescopes shortly after the burst, were reported by Ue- 
mura et al. (2003), Torii et al. (2003), Smith et al. (2003), 
Sato et al. (2003), and Urata et al. (2003). Burenin et al. 
(2003), Bloom et al. (2003) and Matheson et al. (2003) pre- 
sented the results of multi-band follow-up campaigns in the 
optical and near-IR. Greiner et al. (2003b) presented optical 
polarization monitoring of GRB 030329, detecting significant 
variability. These works mostly relied on data collected at 
a single geographical location, thus limiting their ability to 
achieve a continuous temporal coverage. Intensive monitor- 
ing of the afterglow of GRB 030329 in the radio and mil- 
limeter wavelengths was reported by Berger et al. (2003) and 
Sheth et al. (2003), while X-ray observations were reported 
by Tiengo et al. (2003). 

In this paper, we report the results of an intensive, coor- 
dinated, world-wide campaign designed to follow the light 
curve of the OT associated with GRB 030329. Combining 
data from five observatories in three continents, we achieved 
an almost continuous coverage of the OT during the first 
few days after the GRB. Careful cross-calibration was used 
to bring data collected using many different instruments to 
the same reference system. This internally-consistent data 
set has allowed us to correctly incorporate four other sets of 
observations now available in the literature. The final light 
curves which we compiled are unprecedented in their tem- 
poral sampling, and reveal a uniquely rich and complicated 
photometric evolution. We assume throughout a cosmology 
with n,„ = 0.3, D.A = 0.7, //q = 65 km s-^Mpc-'. 

2. OBSERVATIONS AND DATA REDUCTION 

We performed time-resolved CCD photometry of the OT 
at five observatories, from March 29.72 to June 17.19 UT 
(5.66 hr, to 79.7 d after the burst). Observations were car- 
ried out through B, V, R, and / filters, accumulating a total 
of 77, 104, 928, and 96 data points, respectively. Details of 



Table 1. List of observatories 



ID 


Observatory+Instrument 


E. long. 


Filter (number of observations) 


01 


Wise Im+Tek 


34.8 


B(9), V(25), R(484), 1(25) 


02 


Wise Im+Slle 


34.8 


B(8), V(7), R(57) 


m 

UJ 


ooU im+wri 




T3(''^\ A^/">^ T?/^"* Tf^^^ 
ti(J), V(,ZI, KPI, 


05 


MniVT 1 3m4-9 4m4-SITp 


— 1116 




06 


Ptilonitir 1 5ni+Noms 


— 1 16.9 


1(3) 


07 


Mt. Laguna Im+Loral 


-116.4 


B(54), V(69), R(71), 1(67) 


21 


Kyoto" 


^ 135 


CR(391)'' 


31 


RTT'' 


30.3 


B(144), V(167), R(168), 1(165) 


41 


CTIO 1.3" 


-70.8 


B(9), V(13), 1(13) 


51 


FLWO'' 


-110.9 


B(62), V(57), R(lll), 1(57) 


52 


KAIT'^ 


-121.6 


B(14), V(15), R(15), 1(15) 


53 


LCO'' 


-70.7 


B(4), V(4), 1(4) 


54 


LCO-40'' 


-70.7 


R(2) 


55 


KPN04m'' 


-111.6 


B(19), R(4) 


56 


Magellan2'' 


-70.7 


R(25) 


57 


Dupont'' 


-70.7 


B(4) 



Note. — The upper part of the table lists the observatories which took 
part in our campaign. The lower part lists external data sources. 
"Several observatories in Japan. Uemura et al. (2003). 
•"Burenin et al. (2003). 
"Bloom et al. (2003). 
■"Matheson et al. (2003). 

"Unfiltered observations transformed to /?-band, see Uemura et al. (2003). 



Table 2. Reference Stars 



Coordinates (2000) 



Number" V mag" 



Used in sets 



10:44:28.62 +21:27:45.4 


005 


15.587 




2 


3 




6 


10:44:36.86 +21:26:59.1 


016 


13.266 










6 


10:44:39.07 +21:30:59.1 


019 


17.616 


1 




3 




6 ' 


10:44:39.86 +21:34:15.0 


021 


16.842 


1 




3 




6 


10:44:41.75 +21:31:52.6 


026 


19.331 






3 


5 




10:44:42.02 +21:32:32.1 


027 


16.839 




2 


3 


5 




10:44:48.03 +21:34:18.5 


037 


17.909 


1 


2 


3 






10:44:53.66 +21:30:12.1 


047 


18.396 










6 


10:44:54.45 +21:34:29.2 


049 


14.136 






3 






10:44:54.99 +21:29:46.3 




20.275" 








5 




10:44:55.00+21:31:42.9 


050 


19.598 












10:45:06.48 +21:36:13.7 


093 


16.245 




2 








10:45:09.81 +21:35:10.2 


098 


15.469 




2 


3 






10:45:15.36+21:34:16.1 


110 


14.730 




2 


3 




6 



"Star numbers and magnitudes refer to the photometry of Henden (2003). 
Observatory IDs given in TableQ 

"Photometry measured in our dataset, and scaled to the system of Henden 
(2003) 



the equipment and a summary of the observations are given in 
TablelU 

The images were bias-subtracted and flat-field corrected in 
the standard fashion. In each frame we measured the magni- 
tude of the OT, as well as of several reference stars, with an 
aperture of 2 arcsec radius. Three to nine reference stars were 
measured in each subset, depending on the depth and the field 
of view of the images. All of the reference stars were tested 
to be non-variable. The reference stars used in each subset are 
listed in Tabled 

For each of the detectors used, we obtained a set of 
internally-consistent magnitudes of the OT by minimizing the 
scatter of the reference stars over the subset. Outlying mea- 
surements were removed during the process. Apparently un- 
derestimated errors in the photometry of reference stars were 
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increased so that Mof — 1 (where dof is the number of 
degrees of freedom) for each of the reference stars in the sub- 
set. 

Cross calibration, including color terms, of the different in- 
struments was carried out by transforming the photometric 
system of each of the subsets to the one of Henden (2003). 
The transformation parameters were derived by linearly fit- 
ting the differences between the weighted-mean magnitude of 
the reference stars, and their magnitudes in the Henden (2003) 
system, to the colors of the reference stars reported by this au- 
thor Because of the initial uncertainty in the OT colors, two 
iterations were required to transform its magnitude to the ref- 
erence system. The transformation uncertainty of each data 
subset was added in quadrature to the measured photometric 
errors. The photometric errors do not include the uncertainty 
in the zero point of the calibration, 0.02 mag (Henden 2003). 

The photometric measurements obtained at the Mount La- 
guna Observatory failed to yield consistent color terms and 
were therefore cross-calibrated by eliminating the offset be- 
tween overlapping segments of this subset and the rest of the 
light curve. From the residual scatter in these overlapping seg- 
ments, we estimate the resultant systematic cross-calibration 
errors of this set to be ^ 5%. 

We augmented our data set with external photometry from 
Burenin et al. (2003), Uemura et al. (2003), and Bloom et 
al. (2003), kindly provided in digital form by the authors, as 
well as with the data published by Matheson et al. (2003). 
Each of the external data sources was cross-calibrated with 
our set using the procedure described above for the Mount 
Laguna Observatory data. Since the color of the OT changed 
with time, this process may have introduced small systematic 
errors. From the scatter between the external sets and our light 
curves, we estimate this error to be smaller than 5%. 

A few data segments from Uemura et al. (2003), tempo- 
rally overlapping with our own, were not incorporated into 
our light curve because the photometric errors in these seg- 
ments (and hence the scatter of the data points) were greater 
than in our data. We note, however, that our data and those 
of Uemura et al. are in good agreement throughout the over- 
lapping segments. The /-band light curve of Burenin et al. 
(2003), covering the time-span between 0.2 d and 0.6 d, does 
not overlap with our /-band data. The zero point for this seg- 
ment was set by coarsely extrapolating later /-band segments, 
and should therefore be regarded with caution. 

The earliest available observations of the OT (^1.3 hr after 
the burst) were obtained by two groups. A set of unfiltered 
observations were reported by Uemura et al. (2003), who 
transformed the resulting magnitudes to the standard /?-band. 
After cross calibration, these data are in good agreement with 
the rest of our data. A couple of measurements in standard R 
were taken through clouds with the SSO Im telescope (Price 
et al. 2003b). After cross-calibration, the two points yielded 
significantly brighter measurements, by about 0.1 mag. Other 
observations from the SSO Im telescope, taken under better 
conditions, are in good agreement with the rest of the data. 
Because of the good agreement between the bulk of our data 
and those of Uemura et al. (2003), and because of the un- 
favorable weather conditions under which the early SSO im- 
ages were obtained, we preferred the light curve of the for- 
mer over the latter two points. The results we present below, 
however, are insensitive to this selection, and would remain 
essentially the same had we corrected Uemura's data to fit the 
SSO points. The recently reported early observations by Torii 
et al. (2003) are even brighter than the SSO points. Thus, 



Table 3. Observations 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


01 


04 


2729.21443 


1.23024 


16.438 


0.018 


0.050 


01 


04 


2729.21597 


1.23177 


16.414 


0.020 


0.051 


01 


04 


2729.21749 


1.23329 


16.447 


0.018 


0.050 


01 


04 


2729.21900 


1.23481 


16.398 


0.019 


0.050 


01 


04 


2729.22052 


1.23632 


16.415 


0.018 


0.050 


01 


04 


2729.22203 


1.23783 


16.432 


0.018 


0.050 


01 


04 


2729.22355 


1.23936 


16.434 


0.018 


0.050 


01 


03 


2729.22956 


1.24536 


16.786 


0.012 


0.022 


01 


05 


2729.23286 


1.24866 


15.940 


0.011 


0.060 


01 


04 


2729.23549 


1.25129 


16.416 


0.013 


0.048 



Note. — The first ten data points of our light curve. The columns are: (1) 
the observatory code (see TablelT]; (2) the filter used. 2,3,4,5 are for B. V, R, 
I, respectively; (3) JD - 2450000; (4) time since f(,„„, = 2452727.98419757; 
(5) calibrated magnitude; (6) the error in magnitude including the self calibra- 
tion error of the individual observatory; (7) the error in magnitude including 
the uncertainty in cross calibrating the various observatories. 



it appears that forming a consistent picture of the early light 
curve of this event will require further analysis, once all the 
relevant data are available. 

Altogether, we compiled 2687 photometric measurements 
(333, 360, 1644, and 350 in B, V, R, and /, respectively), mak- 
ing GRB 030329 the most extensively studied extragalactic 
cosmic explosion in the optical regime, with the possible ex- 
ception of SN 1987 A. The complete set of photometric mea- 
surements presented in this paper are listed in Table|3l^- 

3. THE LIGHT CURVE 

The final BVRI light curves of GRB 030329 are presented 
in Figure n Our best-sampled set is the /?-band light curve, 
with an almost continuous coverage during 0.05 ^ f < 3 d, and 
a relatively dense sampling up to f « 67 d. We therefore focus 
our attention on this light curve, and, augment our analysis 
with results from other bands when necessary. 

While the OT grossly follows a power-law declining trend, 
it shows remarkably strong short-term deviations from the 
smooth, monotonic decline, particularly between 0.05 ^ t <, 
5 d. Later periods, between 8 < f < 14 d, and at f <; 50, also 
feature pronounced variations of the light curve. We note that 
the light curves do not show a priori any clear distinction be- 
tween "early" and "late" power law components - the signa- 
ture of a break in the decline rate. We also note the lack of a 
clear bump which may be attributed to SN 2003dh, although 
its significant contribution to the total brightness of the OT 
was established spectroscopically. 

At R = 23.25 ± 0.15 mag (Fruchter 2003, Priv. Comm.), 
the host becomes a significant contributor to the observed flux 
during the late decline phase of the OT. Hence below, unless 
stated otherwise, we shall discuss the R-band light curve of 
the OT only, obtained by subtracting the host flux from our 
measurements. The photometric errors of the OT were modi- 
fied to include the uncertainty in the host galaxy magnitude. 

To further discuss the light curve, we decompose it into 
a SN component, and an "afterglow" component, which is 
further divided into two components: a smooth, monotonic 
decline component (the main observational characteristic of 
GRB afterglows), and a perturbations component - variations 

The complete Table with all 2687 measurements is given in the elec- 
tronic version of this paper. The data are also available from our web site: 
http://wise-obs.tau.ac.iVGRB030329/ 
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Fig. 1 . — BVRI light curve of GRB 030329. Our observations are marked by circles, and the external data sets (see TablefTl. cross-calibrated to our photometric 
system, are marked by squares. For presentation purposes the BVI Ught curves are shifted vertically by +0.5, —1.5, and —1.9 mag, respectively. 



over different time scales and intensities about the smooth de- 
cay of GRB 030329. In the following sections we shall con- 
sider each of the OT components separately. 

3.1. The supernova component 

Spectroscopic analysis of SN 2003dh (Stanek et al. 2003; 
Hjorth et al. 2003; Matheson et al. 2003; Chornock et al. 
2003) revealed a remarkable similarity to the well studied SN 
1998bw associated with GRB 980425 (although, see MazzaU 
et al. 2003; Kawabata et al. 2003). In light of this similarity, 
we base our investigation of the SN component in the after- 
glow of GRB 030329 on the known properties of SN 1998bw. 

We constrain the SN component by comparing the OT light 
curve with four alternative models for SN 2003dh. All four 
models are modifications of the light curve of SN 1998bw, 
corrected to the redshift of SN 2003dh (z=0. 1685). The trans- 
formation was carried out by applying synthetic photometry 
(using the methods presented in Poznanski et al. 2002) to 
the large collection of SN 1998bw spectra reported by Patat 
et al. (2001), and taking into account the greater luminos- 
ity distance of SN 2003dh (810 Mpc, compared to 37 Mpc 



of SN 1998bw), cosmological redshift and time-dilation ef- 
fects, as well as the shift in the SN spectrum sampled by 
each filter (K-corrections). The SN light curve was also cor- 
rected for Galactic extinction in the direction of SN 2003dh 
{Eb-v = 0.025 mag, Schlegel, Finkbeiner, & Davis 1998), 
but was not corrected for the host galaxy extinction (see be- 
low). Our derived light curves are in good agreement with 
those computed by Bloom et al. (2003) using a slightly dif- 
ferent approach. 

Comparing the models to the data, we note that any vari- 
ation in the light curve should probably be attributed to the 
afterglow, since no such variations have been detected in the 
light curves of any of the well observed SNe Ic (e.g., SN 
1998bw - Galama et al. 1998b; McKenzie & Schaefer 1999; 
Patat et al. 2001; SN 2002ap - Gal- Yam, Ofek & Shemmer 
2002; Panday et al. 2003; Yoshii et al. 2003; Foley et al. 
2003; SN 1999ex - Sti-itzinger et al. 2002; SN 19941 - Yokoo 
et al. 1994). A valid SN model is therefore required to be 
fainter than the bottom of the dips in the light curve. In partic- 
ular, the minima of two rebrightening episodes in the 7?-band, 
one of ^ 0.3 mag around f « 52 d, and another of ^ 0.2 mag 
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22.5- 



10 20 30 40 50 60 70 80 90 

Time since burst [days] 

Fig. 2. — The host-subtracted S-band hght curve of the OT (open circles) 
compared to four model light curves of SN 2003dh. The plotted lines are: 
the light curve of SN 1998bw, redshifted and K-con'ected (solid line); the 
same model, after applying the best fitting magnitude shift and temporal shift 
relative to the GRB (see text; dotted line); the same model, after applying the 
best fitting magnitude-shift and stretch correction (see text; dashed line); and 
the time-stretched and magnitude-shifted model lowered by 0.3 mag (dash- 
dotted line). 



around f « 68 d (see ^3.3\ set upper limits on the SN mag- 
nitude of 2L93 ± 0.15 and R = 22.04 ± 0. 13 mag at the 
times of these episodes, respectively. 

Figure|2]shows the 7?-band light curve of the OT (open cir- 
cles). Overlaid are the four different model light curves for the 
supernova component. The solid line in Fig.|2]is the redshift- 
corrected light curve of SN 1998bw described above, assum- 
ing it exploded simultaneously with the GRB. This simple 
model is consistently brighter than the data points after day 
~ 20, and is therefore ruled out - had SN 2003dh been iden- 
tical to SN 1998bw, the OT would have been brighter than 
observed after day 20. 

Hjorth et al. (2003) decomposed their observed spectra into 
a SN Ic component (using redshifted versions of SN 1998bw 
spectra) and a power-law component, typical for the optical 
emission from cosmological GRBs. Using this method, these 
authors derived a V-band light curve of the SN component 
during the first month after the GRB. They found that their 
data could not be fitted by a redshifted and K-corrected V- 
band light curve of SN 1998bw (see their Fig. 3), but required 
that the SN component would be slightly brighter (by 0.2 
mag) at the peak, and would also decline much faster than SN 
1998bw, becoming at least 0.7 mag fainter 28 days after the 
GRB. Similar results were obtained by Mazzali et al. (2003), 
who combined spectral analysis with explosion models. 

With this in mind, we consider alternative models for the 
SN component. Following Hjorth et al. (2003) and moti- 
vated by theoretical models of delayed core-collapse (e.g., 
Vietri & Stella 1998; Berezhiani et al. 2003), we consider 
a non-simultaneous SN model. The model light curve (dot- 
ted line) was obtained by fitting the V-band light curve of 
SN 1998bw to that of SN 2003dh (Hjorth et al. 2003), with 
two free parameters: a time lag between the GRB and the SN 
explosion (Ar), and a magnitude-correction (Am). The best 
fit values are: Ar = -4.7+^-^ d and Am 0.08 ±0.10 mag. 



with x^/dof — 2.5/4. This model initially agrees with the 
data, but becomes too bright after day 50, and is therefore 
discarded. Larger temporal offsets (required to make the SN 
fainter at late times) are inconsistent with the spectral analy- 
sis of Stanek et al. (2003), Matheson et al. (2003), and Hjorth 
et al. (2003). In the next model (the dashed line in Fig. |2li, 
we introduced a magnitude correction parameter (Am), and a 
"stretch" parameter s (similar to the formalism used by Perl- 
mutter et al. 1999) which adjusts the width of the light curve 
so that the model luminosity of SN 2003dh at time t after the 
burst is given by the luminosity of the redshifted SN 1998bw 
utt/s. We find s = 0.80 ±0.05 and Am = -0.01 ±0.10 (with 
X^/dof = 0.9/4) by fitting the V-band light curve of SN 
1998bw to the measurements of Hjorth et al. (2003). This 
model provides a somewhat better fit to the data compared to 
the delayed-GRB model (dotted line), but still overestimates 
the data after day 50. 

In search for a SN model which will satisfy all available 
constraints, we consider two additional inputs. Bloom et al. 
(2003) used a self-consistent multicolor data set obtained with 
a single instrument (ANDICAM, mounted on the CTIO 1.3m 
telescope). They found that their data could be adequately 
modeled using the redshifted and K-corrected light curves of 
a SN 1998bw-like event, which is slightly brighter than SN 
1998bw. Such light curves are ruled out by our measurements, 
obtained from day 20 onward. The probable reason for SN 
1998bw being consistent with the observations of Bloom et 
al. (2003), is that their data are well sampled until day 12, 
and their latest data points were obtained at day 23, just when 
our data start indicating that SN 1998bw-like light curves are 
too bright. 

Matheson et al. (2003) did not pursue the full calcula- 
tion of K-corrections required to produce the light curves of 
SN 1998bw as they would appear for an event exploding at 
z = 0.1685. Instead, they used the V-band light curve of SN 
1998bw as a proxy for the /?-band light curve of SN 2003dh. 
Comparing their adopted "7?-band" light curve with the results 
of our own calculation, we find that their light curve has a 
similar shape, but is fainter by ^ 0.4 mag, taking into account 
their estimated extinction Ar < 0.2 mag. In their SN light 
curve model, Matheson et al. (2003) further scaled down SN 
2003dh by 0.2 mag, relative to their adopted SN 1998bw 
light curves. Thus, the good fit they reported for SN 1998bw- 
like light curves is actually for an event fainter by ^ 0.6 mag. 

Because our late-time data rule out non-stretched SN 
1998bw-like light curves, even if they are as faint as advo- 
cated by Matheson et al. (2003), we are led to construct 
a model, which is identical to the stretched model (dashed- 
line), but is attenuated by 0.3 mag (the dash-dotted in Fig.|3 
- a compromise between the results of Bloom et al. (2003) 
and Hjorth et al. (2003), suggesting peak magnitudes at least 
as bright as those of SN 1998bw, and those of Matheson et 
al. (2003), who found that SN 2003dh was fainter than SN 
1998bw by ^ 0.6 mag. This last model seems to best fit both 
our data and other available studies, and is also in good agree- 
ment with the theoretical modeling reported by Mazzali et al. 
(2003). Both the attenuation factor and the stretch parame- 
ter of this last model are well situated in the parameter space 
of fitted 1998bw-like SN, which was obtained by Zeh, Klose 
& Hartman (2003) for a sample of seven GRB optical coun- 
terparts. Note, that we do not correct for the host-galaxy ex- 
tinction of SN 2003dh. If we assume negligible extinction of 
SN 1998bw and the maximal extinction values allowed by the 
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Fig. 3. — R-band light curve of the afterglow. The light curve was derived 
by subtracting the contribution of SN 2003dh and of the host galaxy from 
the observed light curve. Overlaid are a few power-law fits to early and late 
times (see text for details). A break at ? ~ 5 days is apparent. The slope of 
the pre-break power-law is ai 1.1, and that of the post-break is a2 ~ 2.0. 



analysis of Matheson et al. (2003) and Bloom et al. (2003), 
Av ^ 0.4, our adopted peak magnitude may be consistent with 
that of a SN which is intrinsically identical to SN 1998bw. 
The Am = 0.3 magnitude we adopted for the attenuation of 
SN2003dh relative to SN1998bw, is in fact a lower limit on 
the attenuation (or conversely, an upper limit on the allowed 
peak luminosity). Therefore, below we further check the ro- 
bustness of our results to the value of attenuation we adopt in 
our SN model. 

3.2. The afterglow component 

The clearcut connection between GRB 030329 and SN 
2003dh makes it important to characterize the afterglow in 
detail, in order to check whether this is a typical burst, or an 
exceptional one. Figure|3lshows the 7?-band light curve of the 
afterglow, derived by subtracting the host galaxy flux and the 
preferred SN model (dash-dotted line in Figure|2l see 
details). The afterglow light curve is clearly not a smooth de- 
cline - strong variations are apparent, starting very early afte r 
the burst, and continuing throughout our observations ( t|3.3> . 
Inspection by eye shows that the early light curve may be de- 
scribed by a power-law decline, with a steepening of the slope 
around f ^ 5 d. 

We therefore begin our examination by investigating the 
simplest plausible model - a singly-broken power law. We 
fitted sets of double power laws to the afterglow light curve. 
In each fit, we assumed a different break time f^^™™''^ (i.e., 
one power law was fitted to all the points with t < f^rl"™'^'' ^^'^ 
another one to points with t > f^*™™''''). The intersection time 
(fx) of the two power-laws, and the sum of the values of 
the two fits were calculated for each During the first 

three days, the light curve was sampled more intensively than 
at later times. To reduce the overwhelming statistical weight 
of this segment, we diluted it by binning the points. The 
dashed lines in Fig. |3] are a few examples of such fits, with 
'ire™"' = 3,5,8 d. As can be seen, the values of the early- 
and late-time slopes are weakly dependent on the assumed 



break time. We therefore consider these values (ai « 1.1 and 
a2 ~ 2.0, where a is the power-law decay index defined by 
the dependence of flux, /, on time / oc t^"-) to be robustly 
constrained by the data. 

For assumed break times between 1.5 and lid, the cal- 
culated intersection between the two power-law fits fall con- 
sistently between 3 and 5 d. Minimum values were ob- 
tained for assumed breaks between 3 to 6 d. When car- 
ried out with light curves in the B, V and / bands, the same 
procedure yielded a similar range for the intersection points: 
3 <; f X 8 d. 

To test the sensitivity of our results to the SN model which 
was subtracted from the light curve, we repeated the test using 
several different SN models. All the models were redshift- 
and K-corrected light curves of SN 1998bw, stretched by 5 = 
0.80. We varied the values of the attenuation Am between 
0.3 < Am < 0.6 mag, since brighte r val ues (Am < 0.3) are 
ruled out by our light curves (see § 13. 1> . and fainter values 
(Am > 0.6) are inconsistent with the strength of the observed 
SN features in the OT spectra (Matheson et al. 2003; Hjorth et 
al. 2003). The early-time slope did not change significantly, 
while the late-time slope varied between 1.8 — 2.0, leaving the 
range of intersection times unchanged. 

To conclude, for a singly broken power law model, our data 
robustly constrain the early- and late-time decay slopes to lie 
around ai ~ 1.1 and a2 ~ 2.0, respectively. These values 
depend very weakly on the light curve segments used, or on 
the SN model subtracted. If we take the intersection time of 
these slopes as an estimate for the break time, the resulting 
values lie around f ^ 5 d. However, because of the strong 
variations in the light curve, the time of the break is not well 
constrained, and could be placed at any time between 3 to 
^ 8 d. If we adopt t[,reak — 5d for the time of the break, a fit to 
the data before and after the break yields ai = 1 . 1 1 ±0.01 and 
a2 = 1.96 ±0.02. Maintaining self-consistency, the intersec- 
tion of these two power laws is fx 4.9 d. This result agrees 
with the late-time power-law slope, a = 2.05, measured by 
Matheson et al. (2003) from data obtained between f « 5 d to 
f « 61 d. It appears that a double-power law, broken some- 
where between day 3 and ~ 8, and strongly perturbed by a 
series of bumps (Figures |3] and 0} provides a fair description 
of the complex light curve of this OT. We therefore continue 
and analyze the properties of the variable component derived 
from the simple broken power-law model in the next section. 

Numerous works, analyzing short or sparsely-sampled opti- 
cal light curves, have interpreted some of the deviations from 
a smooth decline observed in the OT light curves as manifes- 
tations of specific phenomena predicted by popular relativistic 
jet models. In particular, Berger et al. (2003) propose a model 
that attempts to consistently account for observations in the 
radio, mm, optical and X-ray bands. We review the complex 
structure seen in our superior data in the context of previous 
analyses in section i)4.1l below. and confront an updated ver- 
sion of the Berger et al. (2003) model with our optical obser- 
vations in ii4.2l 

3.3. The variable component 

Subtracting our best double power-law decline model ( ^'i.2\ 
from the afterglow light curve, we derive the residual 7? -band 
light curve of the afterglow. Three segments of the residual 
light curve are shown in Figures0]|6land0 

Figurel^shows the residual light curve during the first eight 
days after the burst. Arbitrarily, we consider the variability as 
a series of bumps, although with no model at hand, this is a 
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mere matter of convenience. The bumps in Fig. |3 are marked 
according to the notation introduced by Granot, Nakar, & Pi- 
ran (2003), who identified four bumps in the light curve of the 
OT, compiled from reports in GCN circulars (A, B, C, and D 
in Fig.|4}. The early bump, K , was identified by Uemura et al. 
(2003), and the minor bumps. A' and C' are introduced in this 
work. The validity of the minor bump C' is uncertain, because 
the peak of the bump is only ^ la brighter than the dip prior 
to the bump, and since the points forming its peak all come 
from the same data set (Matheson et al 2003). 

The bumps before day ^ 8, which possibly occurr ed pr ior 
to the change in the afterglow decline rate (Fig.|4] see ij l4.3l be- 
low) are asymmetric in shape, with an incline that is typically 
^ 2 times as steep as the decline (with the possible exception 
of C', which is poorly sampled). 

Bumps A, B, and C share a strikingly similar overall struc- 
ture and time-scale during ~ 0.6 d about their maximum. This 
similarity is demonstrated in Fig. [S] (panels a, b, and c), where 
each of the three bumps is shown, compared to a curve man- 
ifesting the common coarse structure of these bumps. To 
derive this "standard profile" of the bumps, we superposed 
bumps B and C over bump A. Bumps B and C were shifted 
by Ar -L02 d. Am -0.008 mag, and ^T = -L718 d, 
A;7i = 0.055 mag, respectively, to obtain a good match by eye. 
The data were then smoothed using cubic splines. Finally, we 
derived a profile curve using a high-degree polynomial fit. ^"^ 

The three events are consistent with this "standard profile", 
comprising of a fast monotonic rise, and a slower, complex 
decline. Panels d, e, and /in Fig. |5l show a comparison be- 
tween the "standard profile" and bumps K, A', and D, respec- 
tively. The brightening rate of bump N , as well as its time 
scale, are similar to the standard profile. However, its struc- 
ture is significantly different. In particular, in contrast to the 
concave early decline which A, B, and C seem to share, the de- 
clining branch of K has a convex form. Our fragmented data 
of bump D are consistent with the standard profile, but on a 
time scale that is longer by a factor of ^ 2 — 2.5. Finally, the 
rising branch of minor bump A' seems to be similar to the tip 
of standard profile, but its decline is slower, perhaps because 
of the rising of bump B. 

Another apparent phase of variability, between 8 d and 13 d, 
features three consecutive low-amplitude (A/? « 0.1 mag) 
bumps, with a time-scale of ~ 1 d (Fig. |6|l. However, the 
amplitudes of the later two bumps {t ^ 10.7 and t ^ 11.7 d) 
are within the cross calibration systematic uncertainties, and 
their maxima data all come from a single source (Mount La- 
guna Observatory), and should be therefore treated with cau- 
tion until confirmed by other observations. A broad bump (F) 
spanning 20 d occurred around f « 30 d (Fig.0. The shape 
of this feature is somewhat sensitive to the SN model used. 
Nevertheless, the deviation from the power-law decline dur- 
ing this period persists for any of the SN models which we 
have tested. 

Strong variability is also detectable during the late decline, 
after f ^ 50 d. In particular, a ~ 0.3 mag rebrightening on a 
time scale of 2 d occurred around f w 52 d (the "jitter episode" 
of Matheson et al. 2003; See Fig.|3|l. Our light curve features 
another ^ 0.2 mag rebrightening on a time scale of a few days, 
around f « 64 d. Because the observations tracing this varia- 
tion were all obtained at the same observatory (MDM 2.4m), 
it is unlikely that this is due to some reduction artifact. To con- 

'** A tabulated version of the "stand ard profile" is available in our website, 
|http://wise-obs.tau.ac.il/GRB030329/| 



elude, it appears that successful models of this well-observed 
GRB should ultimately account for strong variations of the 
optical emission, on time scales of hours to weeks, recurring 
over tens of days after the burst. 

3.4. Color evolution 

Figure |8l shows the B — R color evolution of the light curve 
of GRB 030329. The host and the SN light curve were not 
subtracted, in order to reduce propagated errors and keep the 
results model-independent. To derive the B R light curve we 
interpolated the better sampled 7?-band light curve onto the 
times of the B-band light curve. The interpolation uncertainty 
was calculated using the method of Ofek & Maoz (2003). Fi- 
nally, we binned the B — R color light curves in 0.05 d bins. 

The color evolution of this event was studied in detail by 
Matheson et al. (2003) and Bloom et al. (2003). Both groups 
used subsets of the data presented here. Our analysis yields 
consistent results with those previous works. In particular, 
we measure almost constant colors between days 2 — 5, after 
which the emerging SN component (with colors similar to SN 
1998bw) drove the color evolution of the optical emission. 
Furthermore, the earlier coverage of our light curves enables 
us to measure the color evolution of the OT during the first 
day. 

The B R light curve shows a small (^ 0. 1 mag), but signif- 
icant, color variation during the first day after the burst (Fig- 
ure|8l inset). The color of the OT evolved from B — 7? = 0.58± 
0.01 mag at t = 0.28 d (coiTesponding to Pbr = -0.88 ± 0.02) 
to B — R = 0.66 ± 0.02 mag at f = 0.83 d (corresponding to 
Pbk — —1.07 ±0.05), where figji is the spectral energy slope 
defined by, /v °^ (/y is the specific flux). This early evo- 
lution is unlikely to be related to the SN component. To 
date, color evolution has been observed in the light curves 
of GRB 021004, starting 1.5 d after the burst (Bersier 
et al. 2003; Mirabal et al. 2003), and possibly also in 
GRB 000301C (Rhoads & Fruchter 2001). An interpretation 
of the color change of GRB030329 in the context of the rel- 
ativistic synchrotron model, as a manifestation of the cool- 
ing break frequency going through the optical bands (e.g., 
Galama et al. 2003) is discussed below (see ^4-.l\ . 

4. DISCUSSION 

4.1. Early breaks 

The relative proximity of GRB 030329 has prompted spe- 
cial attention by the astronomical community, and several 
works presenting analysis of various data sets have been pub- 
lished so far. Our superior compilation of optical data, as well 
as the privilege afforded by the availability of these earlier 
works, allow us to inspect some previous suggestions in light 
of the newly available data. 

Early analysis of preliminary optical data revealed that 
the slope of the optical decline became steeper around day 
0.5 —0.6 (e.g. Garnavich, Stanek, & Berlind 2003; Burenin 
et al. 2003; Price et al. 2003b). The change in the power- 
law decline index from 1 to ~ 2, seen both in early optical 
data and in sparse X-ray observations reported by Tiengo et 
al. (2003), combined with the achromatic nature of the break 
(Burenin et al. 2003), seemed to support an interpretation of 
this steepening as a "jet break" - the manifestation of a coni- 
cal geometry in the relativistic emitting material (e.g., Rhoads 
1997). 

Figures|3]and|9]show that this simple interpretation does not 
fit the well sampled light curves now available. A model pos- 
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Fig. 4. — The residual light curve, obtained by subtracting our best fit double power-law (see text) from the light curve of the afterglow, during the first eight 
days after the burst. Five strong bumps, along with two possible minor ones, are apparent in the fight curve. 



tulating a steep (a ^ 2) optical decline starting at day ^ 0.5 
would severely under-predict the optical observations from 
day 1.5 onward. Sustaining a model involving a "jet break" 
at day 0.5 requires an additional source of optical emission 
emerging at day ^ 1. Such a model was indeed suggested 
by Granot et al. (2003) and Berger et al. (2003) and will be 
discussed in the next section. 

As can be seen in Figures 0] and |9l as well as in Uemura 
et al. (2003) and Sato et al. (2003), the early optical data 
(i.e., before the suggested break at day ^ 0.5) are not con- 
sistent with a single power-law decline, as predicted before a 
jet break. In particular, the data presented by Uemura et al. 
(2003), Torii et al. (2003), and Sato et al. (2003) require at 
least one additional break to occur around day ~ 0.25. Our 
compilation clearly elucidates this (Fig. |9j. Several authors 
(Torii et al. 2003; Sato et al. 2003) suggested that this addi- 
tional break (with a change in the temporal slope of Aa ^ 0.3) 
may represent the so-called "cooling break", predicted in the 
context of relativistic fireball models (Sari, Piran & Narayan 
1998) to occur as the cooling-break frequency, Vf, passes 
through the optical band (with Aa = 0.25). As seen in Fig.|9j 
the combination of an early cooling break at day 0.25 and a 



jet break at day 0.5 allows a fair representation of the early 
optical data. Thus, successful models of this event should ac- 
count for at least two early breaks. However, it should be 
stressed that the well sampled light curves available before 
0.25 day are not consistent with a smooth power-law decline, 
and show significant wiggles and bumps (we shall further dis- 
cuss this point below). 

Previous studies suggesting an early cooling break relied 
on the analysis of single-band data, and could therefore only 
detect the predicted shift in the temporal decay slope a, as 
discussed above. Our compilation also allows us to inspect 
the color evolution of the OT around the time of the suggested 
cooling break. As reported in 8 )3.41 we detect a shift in the 
optical spectral index of Aj3 = 0. 19 ± 0.05 between 0.28 (our 
earliest available color data) and 0.83 d after the burst. 

Relativistic synchrotron models predict that a cooling break 
sweeping the optical band would manifest itself by a change 
in the color index of A)3 ~ 0.5. If the observed color change 
is indeed due to the cooling break, then considering that ~ 
40% of the expected color change occured between 0.28 to 
0.83 days, and noting that both theory (Sari et al. 1998) and 
observations (Galama et al. 2003) show that the cooling break 
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Fig. 5. — A close up view of the main bumps detected in the early residual light curve, which was obtained by subtracting our best fit double power-law (see 
text) from the light curve of the afterglow.Panels a, h, c, d, e, and/' display bumps A, B,C, S , A', and D, respectively. The light curves in panels d and e were 
binned, to allow convenient browsing. The solid curves are a smoothed fit to the superimposed light curves of bumps A, B, and C, representing the common 
structure of these three bumps (see text for details). The time scale of the curve in panel/was stretched by a factor of 2.3. 
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Fig. 6. — The variations light curve of the afterglow between day 5 and day 
14. 
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Fig. 7. — The variations light curve of the afterglow between day i8day42. 



evolves in time with a power law, we estimate that the cooling break crossed the V band at day 0.25. We note that the 
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Fig. 8. — The binned {B — R) color of the OT during the first 40 days after 
the burst. The vertical axis on the right shows the corresponding spectral 
power-law, jigR- The inset is a blow-up of the first four days. 
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Fig. 9. — A detailed view of the early Hght curve of the OT associated with 
GRB 030329. The light curve distinctly breaks around 0.25 and 0.5 days. The 
data obtained prior to 0.25 days is not well described by a smooth power law, 
and shows significant wiggles. Note that our data rule out further steepening 
of the light curve at day 0.8, with the slope measured during 0.82 — 0.95 days 
consistent with the slope measured from earlier data (0.6 — 0.78 days) and 
certainly below the values a ~ 2.2 reported elsewhere (see text). Around day 
1 the decHne in the optical flux begins to slow, as the rising branch of bump 
A emerges. 



observed steepening of the slope passed through the optical 
band from high frequencies to low ones, and hence demotes 
models predicting an opposite trend. 

Another way to probe the cooling break frequency value 
during later times is through the optical to X-ray spectral 
index, which is predicted to be constant after the cooling 
break moves below the optical bands. We measured this in 
four epochs of XTE and XMM-Newton observations given by 
Tiengo et al. (2003; t = 0.222, 1.24, 37.24, and 60.85 d; the 
first epoch is the weighted mean of four measurements), by 
interpolating the 7? -band afterglow light curves onto the four 
epochs of the X-ray observations. 

The value of slope at the first epoch —0.93 ± 0.02 (at 
f = 0.222), before or during the proposed cooling transition, 
markedly differs from the values obtained at later epochs, af- 
ter the break, -1.04 ±0.05 and -1.06 ±0.13 at 37.24 and 



60.85 days after the burst, respectively. The value measured 
1 .24 days after the burst, — 0.95 ± 0.04, is consistent with both 
early and later values. A coherent picture therefore emerges 
from the optical observations and X-ray data in which the 
color evolution we detected is possibly the result of the cool- 
ing break passage through the optical band around 0.25 days 
after the burst. Naturally, other explanations for the observed 
early color evolution are possible. 

Finally, we consider the report by Smith et al. (2003) 
of a further steepening of the optical light curve (reaching 
a ~ 2.2) around 0.8 day - distinct from, and occurring after 
the proposed 0.5 day break. This could have confirmed the 
prediction by Granot et al. (2003) that the final slope after a 
jet break around 0.5 days should be steeper than the value re- 
ported (a — 1 .9), as expected from the temporal slope change 
due to a jet break, Aa ^ 1, relative to the post-cooHng break 
slope a ~ 1 .2. However, our superior data, including two mu- 
tually consistent data sets from the MDM and Mount Laguna 
Observatory, obtained through standard filters, are not consis- 
tent with the unfiltered data reported by Smith et al. (2003). 
Fig.|9]shows no segment of the early light curve (including the 
period between 0.8 — 1 days) with a power-law slope steeper 
than a = 2. Thus, the theoretically predicted slope of a = 2.2 
was either not reached, or washed out by the emerging bump 
A. 

To conclude, our results support the interpretation of the 
break around 0.25 days, as a cooling break, within the context 
of relativistic synchrotron models. The interpretation of later 
steepening of the light curve (around 0.5 day) as associated 
with a geometric ("jet") break, can be sustained only if an 
additional emission source, dominating the optical flux from 
day 1.5 onward, is invoked. 

It should also be noted that Dado, Dar, & De Rujula (2003) 
also put forward a theoretical prediction for the optical light 
curves of this event, based on their "cannon ball" model. 
Lacking access to the predicted curves, we are unable to con- 
duct a detailed comparison with our data. However, such a 
comparison is certainly warranted, and could be easily con- 
ducted using our publicly available data. 

4.2. The two-jet model revisited 

Berger et al. (2003) analyzed radio observations of 
GRB 030329, and found that the radio data, as well as mm- 
band observations reported by Sheth et al. (2003), could be 
well described by a relativistic jet model. However, the jet 
parameters indicated a "wide" jet, exhibiting the characteris- 
tic jet break about 10 days after the burst. In view of pre- 
vious analysis of early-optical and X-ray data which seemed 
to require a "narrow" jet, breaking around day 0.5, Berger et 
al. (2003) proposed a composite two-jet model, combining a 
narrow ultra-relativistic component responsible for the 7-rays 
and early (f < 1.5 days) optical and X-ray afterglow, with a 
wide, mildly relativistic component responsible for the radio 
and optical afterglow beyond 1.5 days. This model provides 
a good fit to the radio and mm data as well as to the prelim- 
inary optical data used by these authors. The model also de- 
scribes well the early X-ray observations, but under-estimates 
late ones. The wide jet component provides the additional 
source of emission required to account for the optical obser- 
vations after day 1.5, as discussed above. 

The energy derived from the Frail relation (Frail et al 2001; 
Bloom et al. 2003), assuming only a narrow-jet component 
and the updated parameters used by Bloom et al. (2003), falls 
7c7 below the geometrically corrected mean gamma-ray en- 
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Fig. 10. — A compaiison of the updated two-jet model with multi-band 
data. The model and notation are similar to those presented by Berger et 
al. (2003), except for several minor modifications (see text). The various 
symbols denote data points, as labeled, and the soHd curves are the model 
predictions in each respective band. The model R-band Ught curve describes 
the general trends in the observed light curve well, but some discrepancies 
remain - the model becomes too bright around 10 days after the burst, and 
later on becomes two faint. 




10"' 10° 10' 10^ 
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Fig. 1 1. — The residual light curve obtained by subtracting the modified 
two-jet model from the observed i?-band light curve. Note that at f < 1 day 
the model describes the general light curve trends very well, but significant 
wiggles around the smooth trend remain. A phase of strong variations can 
be seen between 1 — 8 days, followed by broader undulation around the pre- 
dicted peak of SN 2003dh, that may be associated with inadequacies in the 
SN model. During later times (> 20 days), the data seem to require yet an- 
other source of optical emission, which is also strongly variable. The inset is 
a magnified view of the strong variations between, day 1 and day 8 



ergy of Bloom et al. (2003). The contribution of the wide-jet 
component introduced by Berger et al. (2003) brings the total 
energy to within 1 a of the mean value of Bloom et al. (2003). 

With our improved optical light curves at hand, we revisit 
the Berger two-jet model. Fi sure [Tol shows a comparison be- 
tween an updated two-jet model with radio data from Berger 
et al. (2003) and X-ray data from Tiengo et al. (2003), as 
well as our R-band light curve. The model is essentially 
the same as the one described in Berger et al. (2003), ex- 
cept for the following modifi catio ns. Firstly, a cooling break 
component at 0.25 day (see ^4-.li was incorporated into the 
model. Additionally, the temporal emergence of the second 
jet at f 1.5 day was set to t (instead of the law used in 
Berger et at. 2003) to account for the abrupt rise which our 
nearly continuous data show (in fact, an even steeper emer- 
gence is probably required). Finally, the 1998bw-based SN 
model used by Berger at al. (2003), which we showed to be 
t oo b right, was replaced by our best SN model as described in 

EU 

The modified two-jet model fairly describes the general 
trends in our well sampled optical light curves. However, 
some discrepancies still remain. Figure [TTI shows the resid- 
uals obtained by subtracting the updated two-jet model from 
the observed /?-band light curve. At early times (< 1.5 day), 
the model fits the trends in the light curve. Nevertheless, sig- 
nificant undulations (with peak-to-peak ^0.1 mag), are de- 
tected about the smooth model (Fig. I12> . Thus, this, or any 
other model based on broken power-law segments, should ul- 
timately be supplemented with a mechanism explaining the 
bumpy nature of the optical emission, starting at the very ear- 
liest times (e.g., even prior to the cooling break). Following, a 
period of strong variations ensues, lasting till day eight. Com- 
paring the inset in Figure ^2 with Figure |4] we see that the 
two-jet model eliminates bump N (which under this interpre- 
tation is the combined effects of the cooling and jet breaks) 
while bump A, associated with the emergence of the second 
jet, is diminished. However, later structure, including bumps 




Time since burst [days] 

Fig . 12. — Same as Fig. II II but zooming in on the period prior to the emer- 
gence of bump A. The residual undulations after subtraction of the Berger et 
al. (2003) two-jet model light curve are clearly evident. 



B, C, and D, remains. Moreover, both the structure of these 
bumps and their peak-to-peak amplitude, are essentially unal- 
tered. Thus, the modified two-jet model does not eliminate the 
need for an additional, strongly variable component of optical 
emission. 

After day eight, a smooth undulation, lasting till day ^ 20, 
is seen in Figure^J Since this is the period in which the op- 
tical flux is dominated by the emission from SN 2003dh, the 
discrepancy seen may indicate that our model light curve of 
SN 2003dh is at fault. This undulation can be almost com- 
pletely removed by making our model for SN 2003dh fainter 
by 0.3 ma g - a t the lower limit of the range we consider plau- 
sible (see ii3.1> . 

Finally, the data clearly require an additional source of op- 
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tical flux in order to explain the late-time (> 20 days) light 
curve. This extra component (in addition to the narrow and 
wide jets, and SN 2003dh) should rise from day 20 and 
through day 79, where our data end, and is also required to be 
strongly variable in order to explain the late time light curve 
"jitters". This extra flux cannot be attributed to the SN without 
making its light curve very different from that of SN 1998bw, 
and highly variable. Further tests for the relative contribution 
of SN 2003dh to the late-time optical flux can be obtained 
from late-time spectroscopy of this event. Interestingly, as can 
be seen in the lower-right corner of Figure [TO] the modified 
two-jet model also under-predicts the late-time X-ray points 
(at t = 37.24, 60.85 day), perhaps indicating a need for an 
extra source of late-time flux in this band also. 

To conclude, as already shown by Berger et al. (2003), the 
two-jet model provides a fair, self consistent description of 
observational data from the radio to the X-ray. However, in 
view of the large volume and complex nature of the multi- 
band data obtained for this event, further investigation of this 
model, perhaps combined with a prescription for the addi- 
tional variations component required, is warranted. 

4.3. Can a single broken power-law provide a good fit to the 

data ? 

We demonstrated above that if we assume a power-law 
decline - the simplest OT model commonly used, the data 
require differ ent slopes for the early- and late-time decline 
slopes (5 13. 2> . The values of the decline indices, a\ ~ 1.1 and 
a2 ~ 2 are robustly constrained by the data. The transition 
between these decline slopes requires a break in the optical 
light curve 3 to ^ 8 days after the burst. Unfortunately, the 
strong undulations superposed on the smooth decline trend 
throughout this period prevent us from determining the accu- 
rate timing of the power-law break. 

It is intriguing to test whether this simple model, which nat- 
urally emerges from the analysis of the optical data, can con- 
sistently provide a reasonable fit to the entire multi-band data 
set gathered for this burst. In this context, we note the follow- 
ing major points. In the optical regime, this simple model for 
the "smooth" evolution of the OT requires additional emis- 
sion component/components to account for the strong flux 
variability detected from hours to weeks after the burst. The 
ubiquity of these undulations suggests that such a component 
cannot be avoided by more complex models which have thus 
far been proposed. Particularly, as shown ( ii4.2> . the two-jet 
model advocated by Berger et al. (2003) accounts for some 
of the more prominent features in the optical light curves 
of this event ( K and A in Fig. @}, but does not account for 
other, equally prominent ones. Thus, without concrete and 
self-consistent models for the mechanism of the light curve 
variability, the amplitude of the undulations about the simple 
broken power-law model for this burst does not seem to argue 
against the validity of the single break model. 

The X-ray coverage of this burst is regrettably sparse. Com- 
paring the available data reported by Tiengo et al. (2003) with 
our optical light curves, and assuming that the X-ray and op- 
tical flux are correlated, as found by Fox et al. (2003) for 
GRB 021004, we expect the optical to X-ray slope in dex to re- 
main approximately constant. As we have shown in ii4.1l this 
is indeed the case, especially if the early X-ray point, obtained 
before the proposed cooling transition, is discarded. We note, 
however, that when over-plotting the X-ray data over the R- 
band light curve, and scaling the X-ray points to match the 
early light curve, the two latest XMM points fall below their 



expected position, by a factor of 6 — 10. The physical signifi- 
cance of this discrepancy is not clear We note, however, that a 
discrepancy of similar magnitude (but of opposite directio n) is 
found also when comparing the modified two-jet model ( ii4.2> 
with the X-ray data (Fig.^| bottom right corner). If any of 
these models is correct, this may suggest a weak evolution in 
the optical to X-ray spectral slope. 

Finally, it appears that the greatest challenge for the simple 
one-break model is accounting for the results of the radio and 
mm observations. As elaborated by Berger et al. (2003), these 
data require the existence of a mildly-relativistic jet, which is 
expected to demonstrate a break in the optical regime around 
day 10 after the burst. Thus, the single break between ^ 3 
to 8 d after the burst implied by the simple model seems 
to be somewhat in odds with the constraints posed by the 
radio and mm data. However, no confidence intervals have 
been determined for the timing of the break, both in the op- 
tical and in radio. In the optical band, the measurement is 
hindered by strong, multiple bumps and wiggles superposed 
on the smooth light curves. The apparent conflict in the de- 
termination of the break time may perhaps be negotiated by 
future analysis, which would also model the variations in the 
optical band (e.g., by fitting a physically motivated model ex- 
plaining the observed optical light curves). We thus conclude 
that the data at hand cannot rule out the empirically motivated 
singly-broken power-law model for the emission associated 
with GRB 030329. Assuming this model with f/jreait ~ 8 day, 
Eisoiy) as found in !|4.2I and an interstellar matter density 
n = 1.8 cm^^ (Berger et al. 2003), the estimated total en- 
ergy of the event, derived using the Frail relation (Frail et 
al. 2001; Bloom et al. 2003), is 4.4 x 10^" erg. This value 
is within 1.5(7 of the updated geometrically corrected mean 
gamma-ray energy value of Bloom et al. (2003). 

4.4. Light Curve Variations 

As shown above, the optical light curves of the OT associ- 
ated with GRB 030329 exhibit strong undulations superposed 
on the overall smooth trends, detectable shortly (hours) after 
the burst, and still apparent many tens of days later. The vari- 
ations seen during the period which is best-sampled by our 
light curves have a typical time scale of ^ 1 day, and show a 
similar asymmetric structure, with a rising branch about twice 
as steep as the declining one. Both the characteristic ampli- 
tude 50%), and the structure of these variations are only 
weakly influenced by our assumptions about the underlying 
smooth behavior. In particular, similar results are obtained 
when we assume either the simple, empirically motivated, 
singly broken power-law model, or the more complex two-jet 
model of Berger et al. (2003; Figures l4land ll II see inset). 

Similar undulations have been previously observed in the 
light curves of OTs associated with GRBs. A short-time vari- 
ation was detected in the light curves of GRB 000301C (e.g., 
Masetti et al. 2000; Rhoads & Fruchter 2001; Berger et al. 
2000; Sagar et al. 2000). The short time scale and achromatic 
nature of the variation led Garnavich, Loeb, & Stanek (2000) 
to suggest this bump was caused by microlensing of the OT 
by a star in a foreground galaxy, while a more mundane origin 
for the bump - nonuniform ambient density - was proposed 
by Berger et al. (2000). The time scale and structure (fast rise, 
slow decline) of the bump detected in GRB 00030 IC, are sim- 
ilar to those seen in our light curves of GRB 030329. This, 
combined with the different structure seen in optical and IR 
light curves of the OT of GRB 000301C (Sagar et al. 2000), 
which also leads to a poor fit of the multi-band microlens- 
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ing model to the observed /?-band light curve (Garnavich et 
al. 2000) suggests that the bump seen in the light curve of 
GRB 000301C is more likely to be intrinsic to the source, 
than the result of the rare cosmic alignment required for mi- 
crolensing. More recently, numerous bumps were detected in 
the Ught curve of the OT associated with GRB 021004 (e.g., 
Bersier et al. 2003; Mirabal et al. 2003; Fox et al. 2003). 
Here too, the temporal structure and time scale of the bumps 
is similar to those seen in our Ught curves of GRB 030329. 

A seemingly exceptional case to note is that of 
GRB 970508, the second burst for which an OT was iden- 
tified. The optical transient of this burst (e.g., Galama et al. 
1998) showed a major rebrightening around one day after the 
burst. Following a short rise lasting less than a day, the optical 
emission underwent a smooth, power- law decline over ~ 100 
days. The relatively sparse sampling of the light curve of this 
burst after day 10 does not allow us to determine the exact na- 
ture of the light curve decline during late phases (e.g., search 
for late "jitter" periods), but the overall smooth structure and 
very large amplitude of this rebrightening appear to be quite 
different from those of the OT associated with GRB 030329 
(or from any other OT so far detected). 

The common explanations for the short time-scale varia- 
tions invoke complex density structures around the burst (e.g. 
Wang & Loeb 2000; Berger et al. 2000; Lazzati et al. 2002), 
inhomogeneous energy disposition within relativistic conical 
blast waves (the "patchy shell" model; Kumar & Piran 2000a) 
and continued injection of energy by the central engine (re- 
freshed shocks; e.g., Rees & Meszaros 1998; Kumar & Pi- 
ran 2000b). Nakar, Piran & Granot (2003) performed a de- 
tailed theoretical investigation of GRB 021004, and showed 
that variants of all these models can fit the data, although they 
preferred a version of the "Patchy Shell" model. 

Granot et al. (2003) next studied the variability of 
GRB 030329, using a preliminary compilation of early opti- 
cal data. Interpreting our detailed light curves in view of their 
analysis, we note the following points. The significant vari- 
ations observed in the light curve before the proposed cool- 
ing break at day 0.25 probably argue against variable density 
models, in agreement with the conclusions of these authors. 
Interpreting the 0.5 day steepening in the light curve as a jet 
break, Granot et al. (2003) ruled out the patchy shell model, 
since it cannot produce strong variations after the entire jet is 
visible (i.e., after the jet break). Since strong variations are 
observed also well after day 10 (e.g., the "jitter episode"), the 
patchy shell model cannot explain the late variations observed 
in the light curve of GRB 030329. However, if the jet-break 
occurred as late as day eight (e.g., as in the single power-law 
model), the patchy shall model may explain the strong early 
variations. 

The refreshed shocks scenario (Kumar & Piran 2002b), pre- 
dicts that if the shocks occurred before the jet spreading, the 
time scale of the variations would be Af ^ f. Granot et al. 
(2003) showed, however, that if the refreshed shocks occurred 
after the jet spreading, the time scale of the variations would 
be A; ~ t^^^. Considering the roughly constant time scale of 
the variations, and assuming a jet-break at day 0.5, Granot et 
al. (2003) favored the refreshed shocks model for the varia- 
tions in the Ught curves of this burst. However, if a jet-break 
occurred after day ~five, the almost constant time scale of 
the variations would be in contrast with the refreshed shocks 
model advocated by Granot et al. (2003). Inverting this ar- 
gument, if refreshed shocks are shown to be the likely mech- 



anism causing the observed variations, a model with a single 
jet break around day five becomes unphysical. 

The similar structure of the undulations seen in the light 
curve of GRB 021004 (Nakar et al. 2003), GRB 000301C 
(Panaitescu et al. 2001) and GRB030329, may hint on a "stan- 
dard" variabiUty mechanism in OTs. Naturally, many more 
cases should be studied in order to confirm this suggestion. 

5. SUMMARY 

We observed the optical afterglow of the nearby GRB 
030329 from five observatories across the globe. We carefully 
cross-calibrated the observations, and augmented them with 
published data. The final compilation of BVRI light curves is 
unprecedented in its temporal sampling, and reveals complex 
structure. 

Decomposing the light curve into host galaxy, SN, and af- 
terglow components, we showed that SN 2003dh, associated 
with GRB 030329, could not have had a Ught curve identical 
to that of SN 1998bw. Instead, the evolution of SN 2003dh 
is better described by making the Ught curve of SN 1998bw 
fainter by 0.3 mag, and with a time scale that is 0.8 times 
shorter We subtracted this SN model from the light curve of 
GRB 030329, and found that the residual light curve is well 
described by a double power-law, with a break point in the 
range of 3 to ~ 8 days. The power-law slope of the Ught 
curve changed from ai « 1.1 to a2 ~ 2.0. These results are 
very weakly dependent on the SN model used. 

The SN, host-galaxy, and power-law subtracted light curves 
of GRB 030329 show strong variations with time scales rang- 
ing from ~ 0.5 hr to ~ 10 d. The early variations (<y 8 
days) which are well covered by our observations, are typi- 
cally asymmetric with their ascending branch about twice as 
fast as the descending branch. Their typical, peak-to-peak 
time scales are 12 — 24 hr. Three of the bumps (A, B, and C) 
are of a similar structure during ~ 0.6 days about their max- 
imum. Later variations are harder to characterize, because of 
the lower frequency of our sampling. Nonetheless, they seem 
to have a longer time scale. Periods of strong variability are 
still evident in the light curve tens of days after the burst. 

We showed that the OT color changed during the first day 
after the burst. In the context of relativistic synchrotron mod- 
els, this supports the suggestions, based on single -band light 
curves, that the "cooling break" frequency passed through the 
optical bands around 0.25 days. 

We discussed previous analysis of this event and found that 
the simple model involving a "jet break" occurring around day 
0.5, proposed by several authors, is unable to account for the 
detailed optical data available for this event. At least two ad- 
ditional emission sources are required in order to sustain this 
interpretation: one to account for the flux observed from day 
one onward, as done by Berger et al. (2003), and another 
mechanism must be invoked to account for the ubiquitous un- 
dulations detected in the light curves throughout the observed 
period, i.e., both before and after the jet break. 

Examining an updated version of the Berger et al. (2003) 
model, we find that it provides a fair description of available 
multi-band data. However, several discrepancies still need to 
be accounted for, probably requiring a self consistent model 
for the mechanism producing the ubiquitous bumps and wig- 
gles in the light curves of this event. A similar effort is re- 
quired in order to test a simpler, empirically-motivated singly- 
broken power-law model we present. 

A comprehensive effort to model the large volume of 
data coUected for this burst, expUcitly accounting for the 
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variability on all time scales, appears like the next chal- 
lenge in making this unique event a key to understanding 
the GRB phenomenon. In this vein, in order to ensure 
the maximal usefulness of our observations to the commu- 
nity, we make all the data available through our web site: 
,http://wise-obs.tau.ac.iyGRB030329/. 
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